Reflection symmetry instability at high spins in i62,i64Yb 



R. G. Nazmitdinov/' ^ J. Kvasil,'^ and A. Tsvetkov'^ 

' Departament de Fisica, Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain 
'^Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 14-1980 Dubna, Russia 
'^Institute of Particle and Nuclear Physics, Charles University, 
V.Holesovickdch 2, CZ-18000 Praha 8, Czech Republic 
(Dated: February 2, 2008) 

A shape evolution of i62,164y)3 jjj yj-^st states is traced using the self-consistent Skyrme Hartree- 
Fock calculations. We found that nonaxial octupole deformations (in particular, Y31 term) become 
favorable at hO, > 0.4MeV in ^''^Yb, while in ^''^Yb a nonaxial quadrupole shape is dominant at 
fast rotation. The cranked Nilsson model and random phase approximation are used to understand 
the dynamics of octupole correlations in both nuclei. We demonstrate that the disappearance of 
one of the octupole vibrational modes in the rotating frame gives rise to the nonaxial octupole 
deformations in ^''^Yb, while the octupole modes are nonzero in ^®*Yb. 
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Shell structure plays a prominent role in the formation 
of ground and excited nuclear states. It is crucially im- 
portant at fast rotation, when nuclear stability against 
fission is determined by a delicate balance of different 
constituents of a nuclear potential. In a geometrical ap- 
proach , l| , where the concept of shape deformation is one 
of the basic cornerstones, effects produced by quadrupole 
degrees of freedom are well understood for various effec- 
tive potentials. Whether octupole degrees of freedom are 
of any importance and how they are manifested,- these 
questions stimulate a noticeable fraction of experimental 
efforts in high spin physics [3| . 

Numerous calculations exploiting the intrinsic reflec- 
tion symmetry breaking of a nuclear mean field predict 
the existence of stable axial octupole deformation in the 
ground state for Ra-Th [Z - 88, iV - 134) and Ba-Sm 
(Z ^ 58, N ^ 88) nuclei, where strong octupole insta- 
bility is due to the coupling of ^15/2 - 39/2, «i3/2 - /7/2, 
hii/2 — d^/2 orbitals [3[. The potential energy surface pro- 
duced by different calculations is, however, quite shallow. 
Experimental nuclear spectra also show a pronounced 
parity splitting at low angular momenta. Rather dynam- 
ical octupole effects (vibrations) are dominant over static 
effects (octupole deformations) at low spins. The increase 
of the angular momentum decreases pairing correlations 
which reduce the octupole interaction, since they couple 
the orbitals with the same parity. With rotation, the 
density of different parity states is enhancing markedly 
near the Fermi level. Indeed, one observes a smooth de- 
crease of the parity splitting with the increase of the ro- 
tational frequency fairly in a few nuclei. These features 
are nicely reproduced within the cranking-|-Hartree-Fock- 
Bogoliubov (HFB) approach with Gogny forces in Ba- 
Sm region [J]. Similar calculations with Skyrme [5] and 
Gogny forces predict a nonaxial Ysi octupole defor- 
mation in light nuclei at high spins. The results based on 
the Skyrme interaction demonstrate the importance of a 
nonaxial ^32 octupole deformation in actinide nuclei at 
fast rotation [3] ■ Notice that rotation induces the contri- 
bution of nonaxial quadrupole and octupole components 



of effective nuclear potentials. 

A wealth of experimental information accumulated in 
the past few years on high spin states in the transitional 
nuclei Z ^ 70, A'^ 90 prompted us to investigate the 
question to what extent dynamical octupole effects ob- 
served at low spins in this mass region might develop 
static octupole deformations at sufficiently high angular 
momenta. This is a main subject of this Letter. 

To analyse the rotational evolution of equilibrium de- 
formations we have employed the HFODD code [8]. In 
the present calculations the pairing correlations are ne- 
glected, since they are not expected to play an important 
role at large rotational frequencies. Our HF calculations 
based on the cranked Skyrme SKP interaction (see de- 
tails in RefF?!) demonstrate that a spontaneous symmetry 
breaking phenomenon occurs in ^^^Yb at large rotational 
frequencies (see Fig.l). The mean field solution with a 
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FIG. 1: (Color online) The rotational evolution of the SkP 
Skyrme mean field energy: pure nonaxial quadrupole (solid 
line) and nonaxial quadrupole-|-octupole (dashed line) solu- 
tions. 

broken refiection symmetry (with nonzero (Q30), {Q31) 
momenta and a small admixture of {Q32) one, see Fig. 2) 
becomes favorable at Ti^ > 0.4MeV, in contrast to the 
solution with a reflection symmetry. In ^^"^Yb both the 
minimal solutions are very close to each other, while the 
pure quadrupole one determines properties of yrast states 
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FIG. 2: (Color online) The rotational behaviour of quadrupole 
{Q2m) and octupole {Qzm) momenta for a broken (solid 
line) and an unbroken (dashed line) reflection symmetry SkP 
Skyrme mean field solutions for ^^^Yb. 



at fast rotation. Hereafter, (...) means the averaging over 
the mean field vacuum (yrast) state at a given rotational 
frequency Vl. We obtain a similar pattern with the SKIII 
parametrisation, which conforms our main finding, how- 
ever, with shghtly higher values of the energy minima. At 
slow rotation 7if2 < O.SMeV (when pairing is important) 
the results should be taken cautiously, while the calcula- 
tions provide a most credible tendency at > 0.3MeV. 

A transparent physical idea that an instability of a nu- 
clear potential with respect to a given deformation im- 
plies a softening of the corresponding vibrational mode 
[i[ enables us to shed light on the mean field results. The 
cranked shell model that incorporates a random phase 
approximation (RPA) represents a powerful tool to illu- 
minate dynamics of various correlations, evolving from 
shape fluctuations (vibrations) to a nuclear shape insta- 
bility. We use a CRPA approach 9] that consists in a 
self-consistent solution of the cranked Nilsson potential 
for the yrast line and analysis of the low-lying excita- 
tions near the yrast line in the RPA. Our consideration 
is based on the Hamiltonian 



Hn =Ho-2^ KNr + V 



(1) 



The term = Hn + -f^add contains the Nilsson Hamil- 
tonian Hn with quadrupole deformations and the ad- 
ditional term that restores the local Galilean invariance 



of the Nilsson potential, broken in the rotating frame. 
To describe different parity states, the interaction V in- 
cludes separable monopole pairing, monopole-monopole, 
quadrupole-quadrupole and spin-spin terms for the posi- 
tive parity (tt = -I-) and dipole-dipole, octupole-octupole 
terms for the negative parity (tt = — ). The chemical 
potentials Ar (t =n or p) are determined so as to give 
correct average particle numbers (Nr). All multipole 
and spin-multipole operators have a good isospin T and 
signature r = ±1 (see the properties of the matrix ele- 
ments in ReflTol). They are expressed in terms of doubly 
stretched coordinates Xi = {uji/ujo) Xi, which ensure the 
self-consistent conditions at the equilibrium deformation. 
All the details of this approach are thoroughly discussed 
in Ref 9. Although the CRPA approach is based on the 
Nilsson potential, it contains static and dynamic effects 
induced by rotation and serves a useful purpose of inde- 
pendent view of the Skyrme results. 

The mode associated with the rotation about the x- 
axis determines the Thouless-Valatin moment of inertia 
Jtv 



Jrv ' 



(2) 



Here, Hn{^=X) ^^'■^ positive signature term of the 
full Hamiltonian; an angle operator $ is the canon- 
ical partner of the angular momentum operator Jx- 
The solution of these equations leads to the definition 
Jtv = detA/detB where the matrix A(B) has a dimen- 
sion n — 10(11) due to a coupling of different operators 
involved in Ha{l=t) i- 

We recall that a comparison 
of the dynamical moment of inertia JZ*-^-* — —(PE/dQ^ 
{E = {Hn)) with the Thouless-Valatin Jrv moment of 
inertia calculated in the RPA provides a faithful test 
of the self-consistency of a cranking microscopic the- 
ory. The equivalence between these momenta of iner- 
tia certainly holds, if one found a self-consistent mean 
field minimum and spurious solutions are separated from 
the physical ones (see exact results for a solvable model 
in Ref 11). The use of the phenomenological rotational 
dependence of the gap parameter near the backbend- 
ing region (see Ref|3), the Is and P terms of the Nils- 
son potential affect this equivalence. However the re- 
sult achieved in our calculations yields a reasonable 
self-consistency (see Fig. 3). We obtain a good agree- 
ment between our results and experimental kinematic 
^il^p,.i^il)) = n^2I/{EM + 1) - - 1)) and the 

dynamic Q^^Jp- ^ = hdl/dfl = 2h/{fl^{I 1)) 
momenta of inertia for the yrast states. Here, E^{I) 
is an experimental energy for the rotational band v 
{y = yrast, (3, 7, ...) [1^. The drastic change of the kine- 
matic momenta of inertia starts slightly earlier and ends 
slightly later than in experiments. It is known that in 
the backbending region the cranking approach is less re- 
liable [l^. Nevertheless, one observes a similar rotational 
dependence of the magnitude of the calculated and ex- 
perimental kinematic momenta of inertia till 0.6MeV 
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in -'^^^Yb, while a good agreement holds till ~ 0.45MeV 
in ^^^Yb (see below). The results for the dynamic mo- 
menta of inertia reproduce with a reasonable accuracy 
the fluctuations in the backbcnding region, in spite of 
the approximations made in our approach. The consis- 
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FIG. 3; (Color online) ^^^yb (left) and ^''''Yb (right). Top 
panels; equilibrium deformations in /3-7 plane as a function of 
the angular momentum I = (Jx) — 1/2 (in units of h) . Middle 
panels: the kinematic Sf'-^'(r2) = (Jx) /Q momenta of inertia. 
Filled square (circle) is used for experimental positive (nega- 
tive) parity lowest states. Different symbols display the exper- 
imental data associated with B1,B2... bands (the band labels 
are taken in accordance with the definitions given in Refill. 
Bottom panels: the dynamic ^'■^'(ri) momenta of inertia are 
compared with the corresponding experimental values (filled 
squares). Experimental values are connected by thin line to 
guide eyes. The calculated values of the kinematic (middle) 
and dynamic (bottom) momenta of inertia for the yrast line 
are connected by a solid line. The dashed line displays the 
behavior of the Thouless-Valatin moment of inertia (bottom) . 
All results are obtained within the cranked Nilsson plus RPA 
(CRPA) approach. 

tency between the mean field and the RPA results was 
achieved by varying the strength constants of the pair- 
ing and multipole interactions in the RPA (see RefQ) to 
fulfill the conservation laws for both signatures and both 
parities 



- 0, 



= 0, 



where = (J^ + iJy) / J 2{Jx) ■ It results in the sep- 



aration of physical modes (excitations) from those that 
are related to the symmetries broken by the mean field. 
Two Goldstone modes are associated with the violation 
of the particle number (for protons and for neutrons). 
The other two Goldstone modes are related to the trans- 
lational and spherical symmetries of the mean field. The 
last equation yields a negative signature solution TiU,, a 
collective rotational mode arising from the symmetries 
broken by the external rotational field. We recall that 
the decoupling of Py and components of the center of 
mass operator from physical modes is uniquely defined in 
the cranking model with a signature and a parity sym- 
metry The constraint satisfaction of the transla- 
tional symmetry related to the operator Px determines, 
therefore, the isoscalar dipole strength solely and also 
guarantees the decoupling of the spurious center of mass 
motion from physical modes. This is especially impor- 
tant for the calculations of i?(i?l)-transitions. For the 
isovector dipole coupling strength we adopt the standard 
value xi = ■KVi/A{f'^) with Vi = 130MeV 0. The self- 
consistent octupole strength constants obtained for the 
anisotropic harmonic oscillator 1^] are used in our cal- 
culations for octupole excitations, taking into account a 
rotational dependence of different mean field values in 
the original expression. By virtue of these constants we 
nicely reproduce the experimental octupole Routhians 
(see below). 

The different response on the rotation can be under- 
stood from the analysis of excited states. To this aim 
we define the experimental excitation energy in the ro- 
tating frame hiOy{Vl)exp = Ryj^) — Ryr{^) as a function 
of the rotational frequency fl Here, the Routhian 

function R,y{il) — E,y{n) — hflI^{Q). The experimental 
states are classified by the parity tt = ± and a quantum 
number a which is equivalent to our signature r (see 
Figs. 4, 6). The positive signature states (r = -1-1) cor- 
respond to a = which characterizes rotational bands 
with even spin in even-even nuclei. The negative sig- 
nature states (r = —1) correspond to a = 1 and are 
associated with odd spin states in even-even nuclei. The 
energy huj^{n)exp can be compared with the RPA results, 
huj,j{Q), calculated at a given rotational frequency. The 
agreement between our results and experimental data in- 
dicates that quadrupole deformation (axial and nonaxial) 
is a major factor that determines the evolution of the 
mean field till HQ 0.6(0.45)MeV in ^^'^Yh (^^^^yb). 

In ^^^Yb the band B4 crosses the band Bl at the ro- 
tational frequency hil — 0.265MeV (a transition point) 
and becomes the yrast one for hil > 0.265MeV. The 
change of the yrast line structure produces the back- 
bending (see Fig. 3). It should be emphasized that the 
potential-energy surfaces are very shallow. The defor- 
mations shown in Fig. 3 should be not taken as precise 
predictions, since even before the transition point the 
yrast states m ay p ossess small 1^ — 3") 7-deformation 
(see also RefslHlll)- 

To elucidate the structure near 
the transition point we trace the rotational evolution 
of the CRPA solutions and the lowest two-quasiparticle 
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FIG. 4: (Color online) ^"^^Yb. The rotational evolution of 
CRPA solutions; the positive (negative) signature ones with 
even (odd) spins are displayed on left (right) panels. Number 
in a circle denotes the RPA solution number : 1 is the first 
f = 1 RPA solution etc. Different symbols display the exper- 
imental data associated with B1,B2... bands. The yrast line 
position denoted as "0" (top, left) is displayed by dashed line 
on all panels. Top panels: positive parity n — + solutions. 
The redundant mode (right) ui^ — fl is denoted as "0" and 
is displayed by dotted line. Bottom panels: negative parity 
n — — solutions. 



poles. At = each quasiparticle orbital (routhian) 
is characterized by asymptotic Nilsson quantum num- 
bers. These numbers are irrelevant, however, in the ro- 
tating case and are used only for convenience. The first 
positive parity and signature RPA solution (see Fig. 4, 
top left panel) may be identified with /3-excitations at 
hfl < 0.2MeV. The increase of the rotational frequency 
leads to a strong mixing between 7- (the second RPA so- 
lution at low spins) and /3-excitations. At the transition 
point from an axial to nonaxial rotations, hil « 0.26MeV, 
a two-quasiparticle neutron component 3/2[651]3/2[651] 
dominates in the quadrupole phonon structure (^ 92%). 
This fact suggests that the alignment of a pair 113/2 is 
the main mechanism that drives the nucleus to triax- 
ial shapes. The first negative signature RPA solution 
(Fig. 4, top right panel) carries the properties of a gamma- 
vibrational mode (with odd spins) till hil « 0.28MeV and 
does not lead to any instability. 

From Fig. 4 (bottom panel, right) one observes that 
the first negative signature and parity RPA (octupole) 
solution tends to zero. Almost a zero energy gap be- 
tween the negative parity band (B5) and the yrast line 



Yb 



LU 

m 

i 



CO 

LU 

m 



0.024 
0.016 
0.008 
0.000 



8000 
4000 




l(-0)i ^lyr 



— l+1(-,1)i^ lyr 
•- 1-1 (-1)1 ^ lyr 



— 1+2 (-,0), ^ I yr 
-- l(-,0)i ^lyr 
•-- 1-2 (-0)1^ lyr 



-| — — I — I — I — I — r 

.0 0.1 0.2 0.3 0.4 



1+3 (-1)1- lyr 

l+1(-,1)i^ lyr 

-- - l-1(-,1)i ^ lyr 

---l^(-,1)i ^ lyr 




~i — <~\ — <~\ — ' — I — r 
0.0 0.1 0.2 0.3 0.4 



m [ MeV ] 



FIG. 5: (Color online) ^^^^Yb. Reduced B{E1)- and B{E3)- 
transition probabilities from the lowest negative parity one- 
phonon states to the yrast states. We used the effective pro- 
ton (neutron) charge e = N/A{—Z/A) for the calculation of 
-B(i?l)-transitions, while only the proton contribution is taken 
into account for the B{E3)- transitions. Left(right) panels 
correspond to the transition with AI even (A/ odd). All 
results are obtained within the CRPA approach. 



is also observed at 'h^l > 0.45MeV. The collectivity of 
the lowest negative parity RPA solutions of both signa- 
tures increases noticeably with the increase of the rota- 
tional frequency. Results for reduced B{E1)- and B{E3)- 
transition probabilities evidently demonstrate the onset 
of octupole correlations at h^l > O.SMeV (see Fig. 5). For 
example, B{E3, AI — 3)-transition probability increases 
from ^ 0.7 W.u. at Hfl « 0.3 MeV to - 4.5 W.u. at 
nn « 0.45 MeV (here W.u is the Weisskopf unit). Sev- 
eral two-quasiparticle components originated from /111/2 
and 37/2 subshells for protons and ?i3/2 and hQ/2 sub- 
shells for neutrons contribute to the collectivity of the 
lowest negative parity one-phonon states. The maximal 
weight of two-quasiparticle components is ~ 65%. All 
these features reveal the nature of a shape transition at 
hil 0.45MeV. The onset of the static octupole deforma- 
tions becomes feasible, since the instability point found 
in the CRPA approach coincides with the result of the 
Skyrme mean field calculations (see Fig.l). These re- 
sults suggest that the octupole deformations are due to 
the octupole phonon condensation at fast rotation. 

In ^^"'Yb (see Fig. 6), there are observed 15 rotational 
bands so far The positive signature family (with 

even spin) is comprised of Bl, B2, B8, B9, BIO, Bll 
bands of positive and B3, B4, B5 bands of negative par- 
ities. The negative signature family (with odd spins) 
includes B2, B12, B13, B14, and B6, B7, B15 of pos- 
itive and negative parity bands, respectively. At hn « 
0.27MeV the band Bll crosses the ground band Bl. This 
crossing gives rise to the backbending (see Fig. 3). Ac- 
cording to our analysis, the first positive parity and signa- 
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FIG. 6: (Color online) ^®^Yb. Similar to Fig.4. 
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FIG. 7: (Color online) ^'^''Yb. Similar to Fig.5. 



ial rotations is determined by two-quasiparticle neutron 
component from 113/2 subshell. A good agreement be- 
tween data and our results implies that quadrupole de- 
grees of freedom dominate in this domain of rotational 
frequency values. The negative parity RPA solutions of 
both signatures (see Fig. 6) decrease with the rotational 
frequency. The RPA results reproduce with a good ac- 
curacy the rotational behaviour of the negative parity 
states, especially, the lowest one. While the octupole 
phonons are of collective nature (the maximal component 
is ~ 75%), the quadrupole deformed field remains stable, 
even at very high spins. The growing octupole collec- 
tivity is exhibited in the increasing dipole and octupole 
transitions from the negative-parity one phonon states at 
nn > 0.3MeV (see Fig.7). In particular, B{E3, AI = 3)- 
transition probability increases from ^ 1.2 W.u. at 
nn « 0.3 MeV to - 4.4 W.u. at hfl « 0.45 MeV. 
In contrast with ^®^Yb, all negative parity solutions are 
nonzero at hil < 0.6MeV. Octupole correlations are un- 
able to break the reflection symmetry in ^^*Yb, which 
elucidates the result of Skyrme calculations (see Fig.l). 

Summarizing, we found that the internal structure 
of the low- lying excited (tt = +, r — +1) band, 
which crosses the ground band, has almost a pure two- 
quasiparticle character originating from 113/2 neutron 
subshell, in agreement with the prediction of the phe- 
nomenological analysis In both nuclei the rota- 

tion produces a profound affect on shell structure induc- 
ing noticeable octupole correlations at high spins. At 
hil > 0.3MeV we predict the formation of octupole bands 
with strong B{E1)- and B{E3)- transitions to the yrast 
line (see Figs. 5, 7). We found that at hfl « 0.45MeV the 
octupole phonon solution vanishes in the rotating frame 
in ^^^Yb. It results in the onset of the degeneracy be- 
tween the lowest negative parity and negative signature 
band and the positive parity and positive signature yrast 
line. This mechanism explains the spontaneous breaking 
of reflection symmetry of the rotating mean field with 
Skyrme interaction (see Fig.l). In contrast, in ^^"'Yb the 
octupole correlations manifest themselves as a low-lying 
octupole vibrations of the quadrupole deformed rotating 
nucleus, in agreement with the Skyrme results. 
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